Classical singularities and Semi-Poisson statistics 
in quantum chaos and disordered systems 
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We investigate a ID disordered Hamiltonian with a non analytical step-like dispersion relation 
whose level statistics is exactly described by Semi-Poisson statistics(SP). It is shown that this result is 
robust, namely, does not depend neither on the microscopic details of the potential nor on a magnetic 
flux but only on the type of non-analyticity. We also argue that a deterministic kicked rotator with 
a non-analytical step- like potential has the same spectral properties. Semi-Poisson statistics (SP), 
typical of pseudo-integrable billiards, has been frequently claimed to describe critical statistics, 
namely, the level statistics of a disordered system at the Anderson transition (AT). However we 
provide convincing evidence they are indeed different: each of them has its origin in a different type 
of classical singularities. 

PACS numbers: 72.15.Rn, 71.30.+h, 05.45.Df, 05.40.-a 



The properties of a quantum particle in a random po- 
tential is one of the most intensively studied problems 
in condensed matter physics since the landmark paper 
by Anderson 0. According to the one parameter scal- 
ing theory, in more than two dimensions, there exists a 
metal insulator transition for a critical amount of disor- 
der. Unfortunately the AT in three and higher dimen- 
sions takes place in a region of strong disorder not di- 
rectly accessible to current analytical techniques. De- 
spite the lack of rigorous analytical results, it is by now 
well established, mainly through numerical simulations, 
that the AT is fully characterized by the level statistics 
and the anomalous scaling of the eigenfunction moments, 
V q = J d d r\ifj(r)\ 2q oc L- D ^ q -^ with respect to the sam- 
ple size L, where D q is a set of exponents describing the 
AT. Eigenfunctions with such a nontrivial (multi) scaling 
are usually dubbed multifractals 2J. Critical statistics 
the level statistics at the AT, is intermediate be- 
tween Wigner-Dyson (WD) and Poisson statistics. Typ- 
ical features include: scale invariant spectrum 0, level 
repulsion and asymptotically linear number variance 

Both level statistics and multifractal properties are 
universal in the sense that parameters such as the slope 
of the number variance or the set of multifractals expo- 
nents D q depend only on the dimension of the system and 
not on boundary conditions, shape of the system, or the 
microscopic details of the disordered potential. However 
the functional form of level correlators as the level spac- 
ing distribution or the number variance may be affected 
by such variables. 

A natural question to ask is whether critical statistics 
can be reproduced by generalized random matrix models 
(RMM). The answer to this question is positive: Critical 
statistics has been found in RMM based on soft confin- 
ing potentials Q , effective eigenvalue distributions 0, H] 
related to the Calogero-Sutherland model at finite tem- 
perature and random banded matrices with power-law 
decay 0. The latter is specially interesting since an AT 
has been analytically established by mapping the prob- 
lem onto a non linear a model. In a certain region of 



parameters, all of above RMM share the same spectral 
kernel, K(s) — (T is a free parameter which 

enters in the definition of the above models). 

In the context of quantum chaos similarities with 
an AT have also been found in a variety of systems: 
Coulomb billiard[Toj. Anisotropic Kepler problem 
and generalized Kicked rotors [12J ■ In all of them the clas- 
sical potential has a singularity and the classical dynam- 
ics is intermediate between integrable and chaotic. In a 
recent letter [T5 | we have put forward a new universality 
class in quantum chaos based precisely on the relation 
between the type of singularity of the classical potential 
and the properties of the quantum eigenstates. Specifi- 
cally, for a certain kind of non-analyticity, it was found 
that the level statistics is described by critical statistics 
and the eigenstates are multifractal as at the AT. 

Similar properties have also been found in pseudointe- 
grable billiards [T3.ll5| . The level statistics of these mod- 
els is very well described by a phenomenological short 
range plasma model whose joint distribution of eigen- 
values is given by the classical Dyson gas with the log- 
arithmic pairwise interaction restricted to a finite num- 
ber k of nearest neighbors. Explicit expressions for the 
level statistics, usually referred to as SP statistics, are 
available for general k. For k — 2, i?2(s) = 1 — e _4s , 
P(s) = 4se" 2s and Y?{L) = L/2 + (1 - e- 4L )/8 where 
R2(s) = —K(s) 2 is the two level correlation function, 
P{s) is the level spacing distribution and S 2 (L) is the 
number variance. Thus SP statistics reproduces typical 
characteristics of critical statistics as level repulsion, lin- 
ear number variance, with a slope depending on k. 

The aim of this paper is to incorporate the systems 
described by SP into the classification introduced in 
which relate quantum properties with singularities of the 
classical potential. We shall see that differences between 
SP and critical statistical stem from the fact that both 
come from different classical singularities. In order to 
proceed we propose a non-analytic Hamiltonian whose 
level statistics is exactly described by semi-Poisson statis- 
tics. It is then shown that differences between critical 
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and semi-Poisson statistics are due to the fact that both 
come from different type of singularity in the classical 
potential. 




FIG. 1: Level spacing distribution for 7 = n/2 and A — 10 
and different system sizes. The agreement with SP is impres- 
sive even in the tail of P(s) (see Inset). 

The models-We start our investigation with a general- 
ized kicked rotor with a a step-like singularity, 

n = ^ + V(q)J25(t~nT) (1) 

n 

with V(q) a step function, V(q) = a; q G [—/?,/?) and 
otherwise V(q) — ; /3 <E [— 7T, n[ and a is a real number. 
By using the method introduced in 0] one can map the 
evolution matrix associated to Eq. onto the following 
ID Anderson model with long-range hopping. 

Hil> i = e i i> i + J2 F (i-3)i>j ( 2 ) 
3 

where ej ~ tan(m 2 ) are pseudo random numbers (for a > 
1) and HH F(m-n) = fjj d(9tan(V'(<9))e- l9 ( m -™) = 
^ sui7(m-n) w ^ n 4 rea j positive constants related to 
(3, a. The case 7 = n/2 describes the high energy limit 
of the interaction of a charged particle with quenched 
Coulomb scatterers of alternate sign. 

In this paper we shall mainly investigate the Hamil- 
tonian Eq. J2J assuming that ej is a random number 
extracted from a box distribution [-W/2, VP/2]. We pro- 
ceed so in order to make an accurate analysis of the 
level statistics necessary for a precise comparison with 
SP statistics. However a detailed study of the Hamil- 
tonian Eq. including the details of the multifractal 
spectrum and the quantum wave-packet evolution, will 
be published elsewhere 

Results- 

We first state our main results: 
1. For A ^> W and any 7 in Eq. [3 the spectrum is scale 
invariant and the level statistics is exactly described by 



semi-Poisson statistics (these findings are in agreement 
with a unpublished calculations of Charles Schmit) . 2. 
There is a transition from SP to Poisson statistics as 
the strength of disorder is increased. 3. In the region 
A W the eigenfunctions arc multifractal but with 
a multifractal spectrum clearly different from the one 
observed at an AT. 

We start by providing analytical evidence that the level 
statistics of Eq. (J2J is described by SP statistics. We 
express the Hamiltonian (JIJ in Fourier space as, 

H = E k i; k +J2Mk,k'W k (3) 

k^k' 

where E k = £r s -^e lkr and A(k,k') 
~k J2 n e«e -m ( fe ~ fe '). We fix 7 = 7r/2 (our findings 
do not depend on 7), after a simple calculation we 
found that E k is not a smooth function (this step-like 
singularity is indeed the seed for the appearance of SP), 
E k = An/2 for k < n and E k = -An/2 for k > n. There 
are thus only two possible values of the energy separated 
by a gap 8 — An. Upon adding a weak (A S> W) 
disordered potential this degeneracy is lifted and the 
spectrum is composed of two separate bands of size ~ W 
around each of the bare points — An/2, An/2. Since 
the Hamiltonian is invariant under the transformation 
A — * —A, the spectrum must also posses that symmetry. 
That means that, to leading order in A (neglecting 1/A 
corrections), the number of independent eigenvalues of 
Eq. is n/2 instead of n. 

We now show how this degeneracy affects the roots 
(eigenvalues) of the characteristic polynomial P(t) = 
det(H — tl). Let Pdis{t) = a-o + ait + . . . a n t n be the 
characteristic polynomial associated with the disordered 
part of the Hamiltonian. We remark that despite of its 
complicated, its roots, by definition, are random numbers 
with a box distribution [— W, W). On the other hand, in 
the clean case P c ; ea „(t) = (t- A) n / 2 {t + A) n /' 2 (n factors 
are not considered). Due to the A — > —A symmetry the 
full (Eq. |2| case Pf u u corresponds with Pd is but replac- 
ing t k factors by a combination (i — A) kl (t + A) k2 with 
k\ + &2 = k. The roots of Pf u u will be in general com- 
plicated functions of A. However in the limit of interest, 
A ^> W — > 00, an analytical evaluation is possible. By 
setting t = ti — A we look for roots t\ of order the unity in 
the A band. We next perform an expansion of the char- 
acteristics polynomial Pf v u to leading order in A. Thus 
we keep terms A n / 2 and neglect lower powers in A. The 
resulting P fuU is given by, P fuU = t™ /2 + a n - 2 ti n/2 ~ 2 /3 + 
an-3ii" /2 ~ 3 / 4 + • • • 2a n/2+1 ti/n + 2a n/2 /(n + 2) where 
the coefficients a n are the same than those of Pdia above 
but only n/2 of them appear in the full case. The eigen- 
values e ■ of Eq. © around the A band are e ■ = A + j3, 
with (3i a root of Pf u n- The effect of the long range 
interaction is just to remove all the terms with coeffi- 
cients ao to a n /2 from the characteristic polynomial of 
the diagonal disordered case. The spectrum is thus that 
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of a pure diagonal disorder where half of the eigenvalues 
have been removed. The remaining eigenvalues are still 
symmetrically distributed (the ones with largest modu- 
lus are well approximated by t max = ^p^) around 

A. That means, by symmetry considerations, that the 
removed ones must be cither the odd or the even ones. 
This is precisely the definition of semi-Poisson statistics. 
In conclusion, the power-law random banded reproduces 
exactly the mechanism which is utilized in the very defi- 
nition of semi-Poisson statistics. We finally mention that 
the only effect of the coefficients 3,4. . . ,n + 2/2 is to 
renormalizc the effective size of the spectrum: ~ 2W for 
diagonal disorder and ~ 2W/y/3 for the Eq. J2J. 




FIG. 2: Number variance for different sizes N, A and 7. Pro- 
vided that A 1 the number variance is given by SP statis- 
tics for any N, A and 7. In the inset P(s) is shown for A — 10, 
N = 800 and different 7. As shown P(s) is not sensitive to 
the specific value of 7. 

The above analytical arguments have been fully cor- 
roborated by numerical calculations. By using standard 
diagonalization techniques we have obtained the eigen- 
values of the Hamiltonian Eq. for different volumes 
ranging from N = 500 to N = 8400. The number of dif- 
ferent realizations of disorder is chosen such that for each 
N the total number of eigenvalues be at least 5 x 10 5 , in 
all cases W = 1. Eigenvalues close to the band edges 
(around 20%) were discarded from the statistical analy- 
sis. The eigenvalues thus obtained were unfolded (by us- 
ing the splines method) with respect to the mean spectral 
density. We first investigate the level statistics in the re- 
gion A^S> W where, according to the analytical findings 
above, semi-Poisson statistics hold. As shown in Fig. 1, 
the level spacing distribution P(s) (including the tail in 
the inset) does not depend on the system size for vol- 
umes ranging from N = 500 to N — 8400. Moreover 
though level repulsion P(s) « s s <C 1 is still present, 
the asymptotic decay of P(s) (see inset) is exponential as 
for a insulator. All these features arc spectral signatures 
of an AT. 

The study of long range correlators as the number vari- 
ance (see Fig. 2), £ 2 (L) = (L) 2 - (L 2 ) = / Q £ p(x)dx - 

Jq (L — s)R,2(s)ds, (p(x) is the spectral density) further 
confirms this point. It does not depend on the system size 



and its asymptotic behavior is linear S 2 (L) ~ 0.51/ L ^> 
1 as at the AT. 

We now compare the level statistics of the Hamiltonian 
Eq. J5J with SP. As shown in Fig. 1 and Fig 2 (inset) we 
could not detect any perceptible deviation from the SP 
P(s) — 4se~ 2s prediction for different 7's and A ^> W. 
The agreement is impressive even for tail of P(s) (inset). 
Also long-range correlators as the number variance follow 
the semi-Poisson prediction ^(L) = L/2 + (1 — e~ 4i )/8 
for different parameter values (see Fig. 2). Deviation for 
small volumes are well known finite size effects. 

A remark is in order, though the above analysis clearly 
show that the level statistics of our model is described 
by SP and share generic features of an disordered con- 
ductor at the AT there are still important quantitative 
differences. Level statistics at the AT depends on the di- 
mension of the space. For instance in 3D (4£>) the slope 
of the number variance is 0.27 (0.41), by contrast semi- 
Poisson predicts 0.5. Clear differences are also observed 
in short range correlators as P(s) [rH fl8| where it has 
been found that, despite their similarities, critical and 
SP have different functional forms. The reason for such 
discrepancy is as follows: 

As mentioned previously, level statistics at the 3D AT 
(critical statistics) is very accurately described by RMM 
Q whose joint distribution of eigenvalues can be consid- 
ered as an ensemble of free particles at finite temperature 
with a nontrivial statistical interaction. The statistical 
interaction resembles the Vandermonde determinant, and 
the effect of finite temperature is to suppress the corre- 
lations of distant eigenvalues. In SP this suppression is 
abrupt, in contrast to critical statistics, where the effect 
of the temperature is smooth. The reason for the dif- 
ferences between SP and critical statistical is thus due 
to the fact that the interaction among eigenvalues is not 
strictly restricted to nearest neighbors (SP) at the AT. 

We now investigate the eigenvector properties of the 
Hamiltonian Eq. [5] We shall sec that, though the 
they are to some extent multifractals, there exist im- 
portant difference with respect to those at the AT. We 
have studied the scaling of the eigenfunction moments 
P q with respect to the sample size L. For multifrac- 
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FIG. 3: P(s) for 7 = tt/2, N = 800 and different A. A 
transition to Poisson statistics is observed in the limit A <C 1. 
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tal wavefunctions, (P q ) = / d d r\ip{v)\ 2( i cx L~ D ^ q ~^ 
where the bracket stands for ensemble average and D q 
is a set of exponents describing the transition. Below 
we show the multifractal dimensions D q (±10% ) for 
A = 10, 7 = 7r/2 obtained by numerical fitting of (log P2): 
D1.5 - 0.36, £> 2 - 0.30,£> 2 .5 ~ 0.28, D 3 ~ 0.26, L> 4 ~ 
0.24, D 5 ~ 0.22, L> 6 - 0.22. For small g, L> 9 depends 
clearly on q however for larger q the dependence is quite 
weak suggesting that may exist a critical q c such that for 
q > q c D q 6 [0, 1] is a constant. This situation would 
correspond with eigenstates which are truly multifractal 
only up to certain scale. Similar results have been re- 
cently reported for certain triangular billiards |15| though 
in this case it was claimed that D q is constant for any q. 
We remark that at the 3D AT the eigenfunction are truly 
multifractal and consequently D q depends explicitly on 
q for any q. 

We now interpret the special properties of eigenfunc- 
tions and level statistics in the context of the original 
non-random Hamiltonian consisting of a kick rotor with 
a non-analytical step-like potential. From the above ar- 
guments it is clear that the Hamiltonian Eq. JIJ avoids 
dynamical localization typical of a smooth potential due 
exclusively to the step-like singularity of the potential. 
In a recent paper we found that certain types of classi- 
cal singularities induce quantum power-law localization 
of the corresponding eigenvectors. For the case of log 
singularities it was explicitly shown that the eigenvector 
were multifractals and the level statistics was given by 
critical statistics. 

The results of this paper shows that SP can be consid- 
ered as the level statistics associated with chaotic systems 
with classical step-like singularities (in 1 + 1 dimensions) 
or with disordered systems with a step-like dispersion 



relation. We can thus unify different intermediate statis- 
tics ('critical statistics' and 'semi-Poisson statistics') in a 
broader classification based on the universal relation be- 
tween classical singularities and level statistics features. 

Finally we would like to discuss briefly two different 
issues. We have observed (see Fig. 3) that as A becomes 
comparable to W the level statistics shifts slowly toward 
Poisson. The level statistics in this region is still scale 
invariant and even for A <C W small deviations from 
Poisson are not negligible. This finding suggests that the 
model is indeed critical for all values of 7 and A. Another 
issue of interest is the robustness of our results under per- 
turbations. We have added a flux to the Hamiltonian Eq. 
(J2J in order to check whether the breaking of time rever- 
sal invariance has any impact on the level statistics. The 
results are negatives, we have only observed the effect of 
the flux in the s — > limit of the level spacing distribu- 
tion P(s) ~ s 2 instead P(s) ~ s. 

In conclusion we have introduced a new class of 
systems with level statistics described by semi-Poisson 
statistics and multifractal wavefunctions. The appear- 
ance of semi-Poisson statistics has been related to a 
step-like singularities of the classical potential and to 
a singular step-like dispersion relation in a disordered 
system. We have discussed similarities and differences 
with critical statistics and claimed that both are part of 
a larger classification scheme. Finally we have discussed 
the transition to Poisson in our model and the effect of 
a flux on the level statistics. 

I thank E. Bogomolny for suggesting me to look at the 
properties of the Hamiltonian Eq. J5J. I acknowledge 
financial support from the Spanish Ministry of Education 
and Culture. 



[1] P.W. Anderson, Phys. Rev. 109, 1492 (1958). 

[2] H. Aoki, J. Phys. 16C, L205 (1983). 

[3] B.I. Shklovskii, B. Shapiro, B.R. Sears, P. Lambrianides 

and H.B. Shore, Phys. Rev. B 47, (1993) 11487. 
[4] V.E. Kravtsov and K.A. Muttalib, Phys. Rev. Lett. 79, 

1913 (1997). 

[5] B.L. Altsfmler, I.K. Zharekeshev, S.A. Kotochigova and 
B.I. Shklovskii, JETP 67 (1988) 62. 

[6] K.A. Muttalib, Y. Chen, et.al. Phys. Rev. Lett. 71, 
(1993) 471; Y. Chen and K. A. Muttalib, J. Phys.: Con- 
dens. Matter 6 L293 (1994). 

[7] M. Moshe, H. Neuberger and B. Shapiro, Phys. Rev. Lett. 
73, (1994) 1497. 

[8] A. M. Garcia-Garcia and J.J.M. Verbaarschot, Phys. 
Rev. E 67, 046104 (2003). 

[9] F. Evers and A.D. Mirlin, Phys.Rev.Lett. 84, (2000) 
3690;E. Cuevas, M. Ortuno et.al. Phys. Rev. Lett. 88, 
(2002) 016401;A.D. Mirlin, Y.V. Fyodorov, et.al., Phys. 



Rev. E 54, 3221 (1996). 
[10] B.L. Altshuler, L. S. Levitov, Phys. Rep. 288, 487 (1997). 
[11] D. Wintgen and H. Marxer, Phy. Rev. Lett. 60, 971 

(1988). 

[12] B. Hu, B. Li, J. Liu and Y. Gu, Phys. Rev. Lett. 82 
(1999) 4224; J. Liu, W. T. Cheng, C.G. Cheng, Comm. 
Theor. Phys. 33 15 (2000). 

[13] A.M. Garcia-Garcia, J. Wang, Phys. Rev. Lett. 94, 
244102 (2005). 

[14] E. B. Bogomolny, U. Gerland and C. Schmit, Phys. Rev. 

E, 59 (1999) R1315. 
[15] E.B.Bogomolny, C. Schmit, Phys. Rev. Lett. 92, 028944, 

(2004). 

[16] S. Fishman, D. R. Grempel, and R. E. Prange, Phys. 

Rev. Lett. 49, 509 (1982). 
[17] I. Varga and D. Braum, Phys. Rev. B 61, (2000) R11859. 
[18] A.M. Garcia-Garcia and J. Wang. In preparation. 



